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THE COSMIC GAMMA-RAY SPECTRUM 


FROM SECONDARY PARTICLE PRODUCTION IN THE METAGALAXY 


, F. W. Stecker* 

ABSTRACT 

The purpose of this paper is to discuss the form and intensity of the 
spectrum of cosmic gamma-rays resulting from the production and decay of 
neutral pi-mesons produced in metagalactic cosmic-ray p-p collisions. It is 
assumed that intergalactic space contains ionized hydrogen gas at a density of 
10" 5 cm -3 as indicated by recent X-ray observations in the 1.5-8 keV region. 

Using the Friedmann solution to the Einstein field equations of general rela- 
tivity as a description of our expanding universe, a discussion is presented of 
the effects of red-shift and spatial curvature on the generation and distortion of 
the local gamma-ray spectrum from the decay of neutral pi-mesons. Numerical 
calculations are presented for the Einstein-de Sitter solution, which is found to 
be an adequate model for these calculations. Various models are presented to 
represent the possible flux of metagalactic cosmic-rays. In calculating meta- 
galactic gamma -ray spectra, the effect of gamma -ray absorption at large red- 
shifts is taken into account. 

A discussion of the results is given. The results indicate that future 
gamma-ray experiments in the 1 - 100 Me V region may yield valuable information 
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relating to cosmology, cosmogeny, and the metagalactic cosmic-ray flux. In 
particular, the metagalactic gamma-ray spectra predicted tend to peak near 70 
(1 + z ) _1 MeV where z ma , the maximum red-shift at which cosmic rays are 
produced, may correspond to the age of the universe at the epoch of galaxy 


formation. 


THE COSMIC GAMMA-RAY SPECTRUM 


FROM SECONDARY PARTICLE PRODUCTION IN THE METAGALAXY 


F. W. Stecker* 
Goddard Space Flight Center 
Greenbelt, Maryland 


INTRODUCTION 

In recent investigations (Stecker, 1967; Stecker, Tsuruta and Fazio, 1968) 
the author made use of recent accelerator and cosmic-ray data to determine the 
details of the cosmic gamma-ray spectrum from the secondary particles pro- 
duced by cosmic-ray collisions in the galaxy. The purpose of this paper is to 
determine the cosmic gamma-ray spectrum from secondary particles produced 
by cosmic-ray collisions in the metagalaxy. This spectrum will differ from the 
galactic (or local) gamma-ray spectrum because most of the generating collisions 
take place at large distances where we are looking back to a time when the uni- 
verse was more compact and collisions were more frequent. These "early" 
gamma -rays will be of lower energy due to the progressive red-shift of the gen- 
eral cosmic expansion. Although various estimates of the flux of these meta- 
galactic gamma-rays have been made (Ginzburg and Syrovatskii, 1964a, b; Gould 
and Burbridge, 1965; Garmire and Kraushaar, 1965), none of these workers have 
taken cosmological factors into account in order properly to calculate a spectrum. 
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THE COSMOLOGICAL EQUATIONS 


For the purpose of these calculations, we may consider models of the uni- 
verse which are both homogeneous and isotropic on a large scale. Such models 
can, in general, be described by the Robertson-Walker line element 


ds 2 - c 2 dt 2 - d-i 2 - c 2 dt 2 - R 2 (t)du 2 (1) 

The time -separable form of the metric is a reflection of the postulated uni- 
formity such that at every moment of world-time the three-space metric is the 
same at all points and in every direction. Such a three-space is a space of con- 
stant Riemannian curvature which may be positive, zero or negative. These 
three alternatives will be designated by k = +1, 0, -1 respectively. If k = +1, 
the universe is closed and finite; if k = 0 the universe is Euclidean, open and 
infinite; if k = -1 the universe is open, infinite and increasingly divergent in 
time. More precisely, 


du 2 


dr 2 

1 - kr 2 


r 2 (d0 2 + sin 2 9 dqt> 2 ) 


so that du measured along the radial direction is given by 


u 


f-jL= 

J 0 Vl - kr 2 


(sin 


- 1 


r 


= < r 


(sinh 1 r 


for k = +1 

for k = 0 

for k = -1 


( 2 ) 
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Gamma-rays travel along geodesics such that ds = 0. Assuming a gamma-ray 
is emitted at a time t e in an interval At e and received at time t r in an interval 
At r , it can be shown that 


u 



constant 


( 3 ) 


and 


CM, B(t,) 
cAt, - R(t.) 

(For a more detailed discussion of the cosmological relations, the reader is re- 
ferred to the excellent articles by Sandage referred to in the text.) The time 
interval, At e , is therefore dilated so that 


JL »(*.) 1 

' *(*.) (4) 
The gamma-ray is therefore shifted to lower energy by an amount 

A\ R ( t r) 

z = x = ¥|t7) " 1 < 5 > 
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We thus obtain the important relation between red-shift and radius of the uni- 


verse given by 


R ( t r) 

R(t e ) = 1 + *(*•) (6> 

It follows from (6) that in a universe where most of the energy density is in the 
form of matter 


n (0 

W7) = (1 + z) 


(7) 


T y (t e ) 


(i +*) 


( 8 ) 


and 


n y (*e) 

n r( tr ) 


= (l + z) ; 


(9) 


where n(t), T y (t) and n y (t) are the average particle density of matter and 
temperature and photon density of cosmic blackbody radiation in the universe. 

We will hereafter designate local (z = 0) quantities with a subscript zero. 
Let f (E r ) be the gamma-ray spectrum generated by the galactic cosmic-ray 
spectrum, I g (E p ) in traveling a unit particle length (1 cm -1 ) through the inter- 
galactic medium. (This spectrum is the same as the quantity l(E y )y^nL> cal- 
culated by Stecker (1967).) 
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We now assume that some ubiquitous generating mechanism causes cosmic- 
rays to be produced with the same power law throughout the universe as observed 
at the earth, so that the metagalactic cosmic-ray spectrum differs only in ab- 
solute intensity from the galactic cosmic-ray spectrum. 11 follows that the form 
of the cosmic gamma-ray spectrum anywhere in the metagalaxy, when observed 
in the co-moving frame at that point, will be the same as the form of f (E y ). 

We may then write down an expression for the integrated metagalactic gamma- 
ray flux in any direction as 



(l + z(^)) 



( 10 ) 


where the factor, (1 + z), takes into account the reduction in flux due to the time 
dilation factor and e~ T represents absorption of gamma-rays along the line of 
sight, I g is the galactic cosmic-ray flux and i(-t) is the cosmic-ray flux at a 
distance -C. Equation (10) may be put into a much more convenient form by ex- 
pressing it as an integral over z. We then obtain 


■W = j 

Jo 


I(z) 

dz n( z) j — 

g 


f (V z ) 

1 + z 


'( E . 


0 d' l 

dz 


( 11 ) 


Since the energy of a gamma-ray is directly proportional to its frequency, it 
follows that 

f(V z) = f[(l+z)E r ] (12) 
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It also follows from (7) that 


n(z) - n 0 (l + z ) 3 


(13) 


The quantity 


d-l du 

dz ®( z ) dz 


(14) 


depends, in general, both upon the cosmological model involved and the epoch of 
world-time which defines the acceleration (or deceleration) of the expansion. 

In Friedmann -type solutions to the Einstein equations, it is found that the ex- 
pansion of the universe is decelerating. The magnitude of this deceleration is 
usually denoted by the deceleration parameter q. In the usual notation, the Hubble 
expansion parameter, H, and the quantity q are defined by the relations 


R(t) 




R(t) 


(15) 


and 


n = - 

q R(t)H 2 

In a decelerating universe, therefore, q > 0. Solutions to the Einstein equation 
with zero pressure and a cosmological constant of zero may be expressed in 
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parametric form in terms of a development angle, 8 (Sandage 1961b) as 

R = a( 1 - cos 8) , , 

t = ^-(8 -sin 8). (16) 

for k = +1, and 

R = a(cosh 8 - 1 ) , 

t - f" (sinh 8- 8) , (17) 

for k = -1, where a = 4^GpR 3 /3c 2 and pis the density of matter in the uni- 
verse (so that pR 3 = constant). 

For the Euclidean case of k = 0, R(t) can be expressed explicitly in terms 
of t by the relation 


R(t) = (6 ttGpR 3 ) 1/3 t 2/3 

From (15) and (18) it then follows that for k = 0 that 

q = y for all t . 

For k = +1, it follows from (15) from (16) that 


(18) 


(19) 


1 - cos 8 
sin 2 8 


(20) 
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For k = -1, it follows from (15) and (17) that q is given by 


1 ~ cosh 1 9 
sinh 2 6 


( 21 ) 


For 6 «1, corresponding to an early epoch of the expansion, (20) and (21) both 
reduce to the Euclidean case of q = 1/2. The Euclidean (Einstein-de Sitter) 
model is therefore a good approximation to the universe if it has not yet reached 
a highly evolved state. It is also compatable with the most probable values of q 
as discussed by Sandage (1961a, 1962), based on the observed magnitude-red 
shift relation, and with the recent determination by Henry, Fritz, Meekins, 
Friedman and Byram (1968) of a mean metagalactic gas density of the order of 
10' 5 cm' 3 . 

Under the assumption of a Euclidean model, we will now determine the 
cosmological effects on the metagalactic gamma-ray spectrum. It can be shown 
(Sandage 1961b) that 


d/ 

dz 


cH 0 -' 


(l ♦ z) 2 (l + 2q 0 Z ) 1 


/2 


( 22 ) 


where cH 0 1 = 10 28 cm. In the Euclidean case, q 0 = 1/2 and we may take in 
equation (11) 


dl 

dz 


10 28 

( 1 + z ) 5/2 


( 23 ) 
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ABSORPTION OF METAGALACTIC GAMMA-RAYS 

An excellent discussion of the absorption processes affecting cosmic gamma- 
rays has been given by Fazio (1967). The principal absorption process to be 
considered is that of electron-positron pair production through interaction with 
the universal black-body radiation field, i.e., the reaction 

y + y -* e + + e _ (24) 

Detailed calculations of the energy-dependent absorption probability for 
this process have been performed by Gould and Schreder (1967). They have 
shown that for a gamma-ray of energy E, y interacting with a black -body radia- 
tion field of temperature l y 



where 


£ 


(me 2 ) 2 

kT E 

y y 


» 


1 


where a - 1/137 is the fine-structure constant, A = B/mc = 3.86 x 10" 11 cm, 
and k here is Boltzmann's constant. The local black-body temperature has 
been found by Stokes, Partridge and Wilkinson (1967) to be 


T 


o 


n 2. 7°K . 


(26) 
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so that the condition £ » 1 corresponds to the condition 


E 


7 


1 ■ 12 x IQ 6 GeV 
(1+z) 2 


(27) 


(see equation (8) and (12)). 

We will restrict ourselves here to a determination of the gamma-ray 
spectrum below 1 GeV and z <10 3 (as will be discussed later) so that the ap- 
proximation given by equation (25) will be generally valid. Therefore, from 
(23) and (25) , we find 

1 . 12 x 10 6 e" 
exp 

r* 1_ (l + y) 2 E y J 

T (V Z ) = 3.9 x 10«E r -« j dy — (28) 


(See appendix for further discussion.) 

THE METAGALACTIC COSMIC-RAY SPECTRUM 

It now remains only to specify a suitable model for the metagalactic cosmic 
ray flux. We will assume that at some early epoch, corresponding to z > z max 
conditions were unsuitable for the acceleration of cosmic-rays. We will consider 
z max correspond to the epoch of galaxy formation and consider two possible 
models for the origin of a metagalactic cosmic-ray flux. For model I, we will 
assume that the metagalactic flux arises through a constant leakage rate from 
the halos of galaxies from z = z max to z = 0. For model II, we will assume 
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that this flux was created primarily in a burst at the time of galaxy formation. 
Thus, model I and model II correspond to the two extreme cases which may be 
expected.* For z max , we will also consider two extremes. One extreme is 
z max = 10 3 > which corresponds to the earliest epoch when galaxy formation could 
probably occur. At z = 10 3 , the black -body temperature of the universe was of 
the order of 10 3 - 10 4o K, cool enough for ionized hydrogen to combine to form a 
neutral gas. According to Peebles (1965), z = 10 3 also corresponds to the epoch 
when gas clouds may begin to form gravitationally bound systems. 

The other extreme for z max which we may consider corresponds to the 
highest red-shift yet observed for a quasar, vis., 2.2. This is, of course, an 
extreme which is technique-limited rather than being limited by any physical 
criteria, and it is included mainly for purposes of discussion. We will also 
consider various intermediate values for z max of 4, 9, and 10 2 . (Doroshkevich, 
et. al. (1967) suggest that galaxy formation took place at z =10-20 whereas 
Weymann (1967) suggests z = 10 2 .) 

It is important to note here that the upper limit, z max , may be effectively 
restricted, not by the epoch of galaxy formation, but by attenuation of the meta- 
galactic cosmic-ray flux due to the collisions themselves. The cross-section 
for inelastic cosmic-ray p-p collisions is of the order of 30 mb. Therefore, 
the lifetime of the metagalactic cosmic-rays against collisional losses is given 

*Models allowing for the possibility of increased cosmic-ray production by galaxies in the past, as 
well as other possibilities, will be considered in a future paper. In particular, it may be of inter- 
est to examine the possibility that the isotropic X-ray spectrum may be a red-shifted 77 “-gamma- 
ray spectrum, a hypothesis which may explain both the intensity and power-law form of the spectrum. 


11 


by 


r 

c 


1 

ncrc 


1 

n Q cc ( 1 + z) 3 


- 10 20 (1 + z) 3 sec 


( 29 ) 


for n Q = 10 5 cm' 3 . 

The lifetime of the universe at a red-shift z is given by 

r u - 10 17 (1 + z)" 3/2 sec (30) 

Cosmic-rays cannot accumulate in the metagalaxy if the ratio, r c / t u < 1. 
The condition r /r = 1 therefore defines a critical value of z = 10 2 beyond 
which a further buildup of metagalactic cosmic-rays cannot occur. With these 
limitations on z in mind, we will now consider the various ideal models for de- 
scribing the metagalactic cosmic-ray flux. 

Model I: In this model, we assume that the galaxies were created on the 
order of 10 10 years ago, as indicated by stellar evolution studies in our own 
galaxy. We will assume that ours is a typical galaxy for which the leakage 
time of cosmic rays from the halo is on the order of 10 8 years. Therefore each 
galaxy has emitted about 10 2 halo-volumes of cosmic-rays since the time of 
galaxy formation. The ratio of the volume of one galactic halo to the volume of 
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the universe is 



( 31 ) 


Taking the cosmic-ray flux in our galaxy as an average for all galaxies and 
taking the number of galaxies in the universe to be 3 x 10 9 (Allen, 1963), we find 
that 


I 0 * 10“ 16 x 3 x 10 9 x 10 2 I g = 3 x l(T s I g (32) 

We assume a constant leakage rate so that the total number of cosmic -rays 
in the metagalaxy is proportional to the time elapsed since galaxy formation. It 
follows from (18) that this time is given by 

T g - 10 10 [(1+ z)" 3/2 -(l + z max )~ 3 / 2 ] yrs . (33) 

The cosmic-ray density will then increase with red-shift according to the 
relation 


I' (g) 

I 


3 x 10" 5 (1 + z) 3 [(1 + z)- 3/2 - (l + z max ) _3/2 ] 


(34) 


However, the cosmic-rays which produce the neutral pi-mesons necessary for 
gamma -ray production are only those above a threshold energy, E h , of about 
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300 MeV (Stecker, 1966). We must therefore determine 


1(0 = I ' (E > E th ; z) 

For a power low cosmic-ray spectrum of the form 


I(>E) ~ E" 3/2 

it follows from the red-shift relation that 


l(E>E,„i z) = 



so that we must use an effective flux of 


( 35 ) 


(36) 


(37) 



3 x 10" s (1 + z) 3 


1 + z 


1 + z_ 


3/2 




(38) 


Model II: In this model we assume that the metagalactic cosmic-rays were 
created in a burst at the time of galaxy formation. Assuming that the total 
number of cosmic-rays released in the initial burst is equal to the total number 
released over 10 1 0 years in model I, we find 


I( z ) 

I 


3 x 10' s (1 + z) 3 



(39) 
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Model III: For a final comparison, we compute the integrated gamma-ray flux 
generated in the galaxies themselves and determine their contribution to the 
metagalactic -gamma -ray flux. We take the average amount of matter in galaxies 
to be about 1% of the total matter in metagalactic space and assume that on the 
average, a fraction of 5 x 10" 2 of this matter is in the form of gas (Allen, 1963; 
Roberts, 1963). Then taking I(z) - I , we find that in this case 

<I(z) n(z)> 

f“ = 5 x 10” 4 (1 + z) 3 (40) 

g o 

Using the models defined by equations (34), (38), (39) and (40), together with 
equations (11), (12), (13), (23) and (28), we have calculated the metagalactic 
gamma-ray spectra produced by models I, II and III. These fluxes are given in 
Figures 1, 2 and 3 respectively. Figures 1-3 also show the gamma ray flux 
expected from the galactic halo in the direction of the pole, taking <nL> = 3 x 10 20 
cm -2 and based on previous calculations (Stecker, 1967).* It can be seen that the 
local gamma-ray spectrum from the galactic halo can be distinguished from the 
metagalactic gamma-ray spectra because the latter are red-shifted and peak at 
lower energies. The metagalactic gamma-ray spectra tend to peak near 7 x 10" 2 / 
(1 + z max ) GeV, being weighted toward higher red-shifts by the effect of greater 
densities at earlier epochs. Because of the density effect, a cosmic-ray burst 


* Thus W < E r > = 3 x 10 : 20 f(E r ). 
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at large red-shifts is much more effective in producing gamma -rays then a 
continuous production of the same number of cosmic-rays. In Figures 1-3, 
we have also indicated the experimental upper limit on the gamma -ray flux 
from the Explorer XI data (Kraushaar, Clark, Garmire, Helmkin, Higbie and 
Agogino, 1965), assuming an integral flux above 0.1 GeV of 3 x 10“ 4 cm -2 sec 
sr” >. 


CONCLUSIONS 

Present evidence about the flux of cosmic -rays between the galaxies is 
quite meager. The most promising way to study the flux is by a satellite ex- 
periment measuring the isotropic gamma-ray flux in the region between 1 and 
100 MeV. Such gamma-rays can supply us with direct information on meta- 
galactic cosmic-rays, because they travel to us in straight lines and suffer little 
absorption. Theoretical metagalactic gamma -ray fluxes from tt 0 decay are pre- 
sented here under various assumptions as to the metagalactic cosmic-ray flux. 
These predictions indicate that an experimental determination of the isotropic 
gamma-ray spectrum at high galactic latitudes and in the energy range 1 - 100 
MeV, could supply valuable information, not only about metagalactic cosmic-rays, 
but also about such fundamental questions as when the galaxies were formed, 
since the metagalactic gamma-ray spectrum will peak near 70 (1 + z max )' 1 MeV. 
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APPENDIX 


GAMMA-RAY OPACITY IN AN EINSTEIN-DE SITTER UNIVERSE 

It has been shown in the text that the gamma-ray absorption from pair 
production through interaction with the universal radiation field is given by the 
factor e“ T where 


(E y , z) = 3.9 x 10 8 E 


8 v ~ 1/2 


f 

Jo 


dy(l+y) 1/2 exp 


1. 12x IQ 6 

(i + y) 2 E. 


- (A-l) 


with E y expressed in GeV. As an intermediate solution to the problem considered 
in the text, numerical solutions where obtained for the implicit relation 


T(E r , z y ) - 1 (A-2) 

which defines the red-shift, z y , beyond which the universe becomes opaque to 
gamma-rays of local energy E y . It was found that the numerical solution to 
equation (A-2) may be quite well approximated by the relation 


1 + z y ^ 2.30 X 10 2 E y 0 ' 484 . (A-3) 

Since the age of the universe corresponding to a red-shift z y , is given by 


t y - 10 10 (l + z y ) _3/2 years , 


(A-4) 
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the earliest epoch from which gamma-rays of energy E y 
formation is found from (A-3) and (A-4) to be 

t y - 2.9 x 10 6 E y 3/4 years . 


can supply us with in- 


(A-5) 


20 


FIGURE CAPTIONS 


Figure 1: Metagalactic gamma-ray spectra from cosmic-ray p-p interactions 

based on a cosmic-ray flux produced by constant leakage from other 
galaxies (Model I) and shown for various maximum red-shifts as 
discussed in the text. Also shown are the local gamma-ray spectrum 
from the galactic halo in the direction of the pole and the upper limit 
implied by the Explorer XI gamma-ray experiment. 

Figure 2: Metagalactic gamma-ray spectra from cosmic -ray p-p interactions 

based on a cosmic-ray flux produced by a burst of cosmic rays at 
z (Model II) as discussed in the text. Also shown are the local 
gamma -ray spectrum from the galactic halo in the direction of the 
pole and the upper limit implied by the Explorer XI gamma-ray 
experiment. 

Figure 3: Metagalactic gamma-ray spectra from the superposition of gamma- 

ray spectra produced in all the galaxies taking red-shift, density and 
curvature effects into account as explained in the text (Model III) . 
Also shown are the local gamma-ray spectrum from the galactic 
halo in the direction of the pole and the upper limit implied by the 


Explorer XI gamma -ray experiment. 
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